Neurolmage: Clinical 5 (2014) 128-140 



Neurolmage: 

CLINICAL 




Contents lists available at ScienceDirect 

Neurolmage: Clinical 



journal homepage: www.elsevier.com/locate/ynicl 



Prospective longitudinal MRI study of brain volumes and diffusion 
changes during the first year after moderate to severe traumatic 
brain injury 

Veronika Brezova ab , Kent GOran Moen c d , Toril Skandsen c e , Anne Vik cd , James B. Brewer f h , 
0yvind Salvesen g , Asta K. Haberg b c * 

3 Department of Circulation and Medical Imaging, Norwegian University of Science and Technology (NTNU), Trondheim, Norway 

^Department of Medical Imaging, St. Olav's Hospital, Trondheim University Hospital, Trondheim, Norway 

c Department of Neuroscience, Norwegian University of Science and Technology (NTNU), Trondheim, Norway 

d Department of Neurosurgery, St. Olav's Hospital, Trondheim, Norway 

^Department of Physical Medicine and Rehabilitation, St. Olav's Hospital, Trondheim, Norway 

f Department of Radiology, University of California San Diego, San Diego, USA 

^Department of Cancer Research and Molecular Medicine, Norwegian University of Science and Technology (NTNU), Trondheim, Norway 
h DepartmentofNeurosciences, University of California San Diego, San Diego, USA 



ARTICLE INFO 



ABSTRACT 



Article history: 

Received 14 January 2014 

Received in revised form 14 March 2014 

Accepted 22 March 2014 

Available online 28 March 2014 



Keywords: 

Post-traumatic amnesia 
Diffuse axonal injury 
Glasgow Coma Scale 
ADC 

Outcome 



The objectives of this prospective study in 62 moderate-severe TBI patients were to investigate volume change in 
cortical gray matter (GM), hippocampus, lenticular nucleus, lobar white matter (WM), brainstem and ventricles 
using a within subject design and repeated MRI in the early phase (1-26 days) and 3 and 12 months postinjury 
and to assess changes in GM apparent diffusion coefficient (ADC) in normal appearing tissue in the cortex, hip- 
pocampus and brainstem. The impact of Glasgow Coma Scale (GCS) score at admission, duration of post- 
traumatic amnesia (PTA), and diffusion axonal injury (DAI) grade on brain volumes and ADC values over time 
was assessed. Lastly, we determined if MRI-derived brain volumes from the 3-month scans provided additional, 
significant predictive value to 12-month outcome classified with the Glasgow Outcome Scale— Extended after 
adjusting for GCS, PTA and age. 

Cortical GM loss was rapid, largely finished by 3 months, but the volume reduction was unrelated to GCS score, 
PTA, or presence of DAI. However, cortical GM volume at 3 months was a significant independent predictor of 12- 
month outcome. Volume loss in the hippocampus and lenticular nucleus was protracted and statistically signif- 
icant first at 12 months. Slopes of volume reduction over time for the cortical and subcortical GGM were signif- 
icantly different. Hippocampal volume loss was most pronounced and rapid in individuals with PTA > 2 weeks. 
The 3-month volumes of the hippocampus and lentiform nucleus were the best independent predictors of 12- 
month outcome after adjusting for GCS, PTA and age. In the brainstem, volume loss was significant at both 3 
and 12 months. Brainstem volume reduction was associated with lower GCS score and the presence of DAI. 
Lobar WM volume was significantly decreased first after 12 months. Surprisingly DAI grade had no impact on 
lobar WM volume. Ventricular dilation developed predominantly during the first 3 months, and was strongly as- 
sociated with volume changes in the brainstem and cortical GM, but not lobar WM volume. 
Higher ADC values were detected in the cortex in individuals with severe TBI, DAI and PTA > 2 weeks, from 
3 months. There were no associations between ADC values and brain volumes, and ADC values did not predict 
outcome. 

© 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 



* Corresponding author at: Department of Neuroscience, Faculty of Medicine, 
Norwegian University of Science and Technology (NTNU), Post Box 8905, Trondheim 
7491, Norway. 

E-mail address: asta.haberg@ntnu.no (A.K. Haberg). 



General cerebral atrophy is a common consequence of moderate to 
severe traumatic brain injury (TBI) (Trivedi et al, 2007; Ding et al, 
2008; Sidaros et al, 2009; Warner et al, 2010a). This general atrophy 
is not predicted by focal lesion volume (Marcoux et al, 2008), and 
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develops over at least 3 years with the bulk loss occurring in the first 
year (Sidaros et al, 2009; Blatter et al, 1997). However, there is limited 
knowledge on the trajectories of volume changes for different brain 
structures, especially in humans. There is evidence from rodent models 
of closed head injury that brain region and injury type determine the 
histopathological response, and hence degree of tissue loss and out- 
come (Fox et al., 1998; Immonen et al., 2009; Meaney and Smith, 
2011). Longitudinal studies of structural changes following TBI in 
humans have so far only contained two, often poorly defined, time 
points of MRI (Ross, 201 1 ), which is insufficient to describe trajectories 
of volume loss in different brain structures. Furthermore, the impact of 
different clinical characteristics related to TBI subtype classification, e.g. 
injury severity and presence of diffuse axonal injury (DAI) /traumatic 
axonal injury, on trajectory of volume changes in different brain struc- 
tures has not been explored systematically. 

Based on the animal literature we predicted that cortical volume de- 
clines first, mainly between the early and 3-month scans. To the best of 
our knowledge there are no human studies documenting longitudinal 
cortical volume changes postinjury. Rodent TBI studies consistently 
show that the histopathological response including neuronal necrosis 
and apoptosis is largely completed within 3 months after injury in the 
cortex (Conti et al., 1998; Rodriguez-Paez et al, 2005; Sato et al, 
2001), although volume loss may continue for 12 months (Smith 
et al., 1997). 

In the hippocampus, on the other hand, protracted neuronal loss 
lasting for at least 12 months is described in animals (Conti et al., 
1998; Smith et al., 1997; Williams et al., 2001). In addition, prominent 
changes in synapse morphology and reduced hippocampal 
neurogenesis are reported (Gao et al, 2011; Scheff et al, 2005; Gao 
et al., 2008). Hippocampal volume loss is reported to be greatest during 
the first weeks, but continues for at least 12 months in rodents 
(Immonen et al, 2009; Smith et al., 1997). Thus, we predicted that hip- 
pocampal volume loss in humans following moderate-severe TBI would 
be protracted, and continue during the entire observation period. 

For the lenticular nucleus less information is available, but quite dis- 
parate histopathological results are reported for other subcortical GM 
nuclei. In animals and humans, both very early and slowly evolving his- 
topathological changes have been shown (Hicks et al., 1996; Greer et al., 
201 1 ; Tong et al., 2002). Patient studies have documented both volume 
loss (Warner et al., 2010b; Bendlin et al, 2008; Wilde et al., 2007) and 
no volume change in the lenticular nucleus (Salmond et al, 2005). 
Taken together the data suggest that lenticular nucleus volume most 
likely declines slowly, and we expected to see a significant reduction 
at 12 months. 

Lobar WM volume changes are described as the slowest and take 
place over a very protracted time period due to Wallerian/Wallerian- 
like degradation being lengthy and incomplete in the brain com- 
pared to the peripheral nervous system (Vargas and Barres, 2007; 
Coleman, 2005; Coleman and Freeman, 2010). Significant lobar 
WM loss was therefore expected to develop between 3 and 
12 months after injury. 

Human studies have shown that the brainstem volume is reduced in 
the chronic phase after TBI (Sidaros et al., 2009; Kim et al., 2008), but an- 
imal studies have not focused on this region. We anticipated that 
brainstem volumes decreased slowly, similar to the lobar WM. 

TBI leads to ventricular dilation in both humans and animals. One 
patient study reports that the period of ventricular dilation lasts until 
2-7 months postinjury with no significant increase thereafter (Blatter 
et al, 1997). In rodents ventricular enlargement is described to take 
place until 6 months postinjury (Bouilleret et al., 2009; Liu et al, 
2010). We therefore expected the ventricular volume increase to take 
place mainly between the early and 3-month scans. Ventricular dilation 
is considered to result primarily from WM volume loss (Anderson and 
Bigler, 1995), but if ventricular dilation develops early after TBI, GM vol- 
ume loss may also contribute significantly. We therefore investigated 
the relationship between ventricular volume change over time and 



the volume changes in cortical GM, lobar WM and brainstem to estab- 
lish the relative importance of volume loss in these structures to ven- 
tricular dilation. 

Furthermore, the impact of clinical measures of TBI mechanism 
and severity, i.e. Glasgow Coma Scale (GCS), duration of post- 
traumatic amnesia (PTA), and DAI, on volume loss in different 
brain structures was assessed. Lower GCS scores were expected to 
lead to more severe loss of lobar WM and brainstem volumes, as un- 
consciousness postinjury is associated with deeper lesions (Jenkins 
et al., 1986) and increasing mechanical forces converging in the mid- 
brain (Meaney and Smith, 201 1 ). Longer duration of PTA (>2 weeks) 
was predicted to be associated with lower hippocampal volumes 
since PTA is a manifestation of hippocampal dysfunction and hence 
likely degree of injury. DAI was predicted to result in increased loss 
of lobar WM volumes, and increasing DAI severity was predicted to 
lead to greater brainstem volume loss. Based on postmortem evi- 
dence, DAI also gives rise to increased cortical thinning and could 
thus reduce cortical volume (Maxwell et al, 2010). 

GM volume changes could be accompanied by changes in GM micro- 
structure. The apparent diffusion coefficient (ADC) obtained from diffu- 
sion weighted imaging (DWI) is a measure of the degree of random 
isotropic diffusion. In GM diffusion is not restricted to a certain direction 
as in WM, and ADC values are considered to reflect cell number (neu- 
rons and for instance inflammatory cells), cell size (swelling vs. shrink- 
age), as well as other cytoarchitectonic features such as the composition 
and fraction of neuropil and extracellular matrix. It has previously been 
demonstrated that TBI patients with lesions visible on conventional MRI 
have increased ADC values in the brainstem the first month after TBI 
(Brandstack et al, 2011), while ADC values in cortex and subcortical 
GM increased between ~2 and ~3 months in moderate-severe TBI 
(Bendlin et al., 2008). Thus we expected the severe TBI patients to 
have increased ADC values in the normal appearing midbrain at 
3 months. On the other hand, the ADC values in normal appearing cor- 
tical GM should increase more between 3 and 12 months in both mod- 
erate and severe TBI. Since DAI leads to loss of neurons, neuronal soma 
shrinkage, axonal changes and myelin loss in cortex (Greer et al, 2011 ; 
Coleman and Freeman, 2010; Maxwell et al., 2010), we expected that 
cortical ADC values would be higher in DAI than in non-DAI patients 
after 3 and 12 months. Moreover, we predicted that longer PTA would 
be associated with increased hippocampal ADC values due to increased 
cell loss. 

The primary objectives of this prospective study in moderate and se- 
vere TBI patients were to investigate the trajectories of volume change 
in total cortical GM, hippocampus, lenticular nucleus, lobar WM, 
brainstem and ventricles using a within subject design and repeated 
MRI in the early phase and 3 and 12 months after injury and to relate 
the volume changes to clinical measures of TBI mechanism and severity; 
GCS score, PTA, and DAI. Changes in local diffusion properties in GM 
were studied to identify any associations between changes in tissue mi- 
crostructure as described with ADC values and the clinical measures of 
TBI mechanism and severity; GCS score, PTA, and DAI. In addition the 
relationship between ADC values and brain volumes was examined. 
Our final aim was to determine if MRI-derived brain volumes and ADC 
values obtained after 3 months provided additional, significant predic- 
tive value to 12-month outcome after adjusting for GCS score, PTA and 
age. 



2. Materials and methods 

The study was approved by the Regional Committee for Medical Re- 
search Ethics and the Norwegian Social Science Data Services. Written 
consent was obtained from the patient or, for individuals underaged 
or incapacitated, their next of kin. Permission was obtained from the 
Norwegian Directorate of Health to use data from the deceased without 
consent from their next of kin. 
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2.1. Subjects 

124 moderate to severe (GCS score < 13) TBI patients aged 
14-65 years were admitted to the Neurosurgical Department, St. Olav's 
Hospital, Trondheim University Hospital, Norway, in the period October 
10, 2004-October 10, 2007. St. Olav's Hospital is the only level 1 trauma 
center in a region of 680,000 inhabitants. 1 7 of these patients did not 
participate in MRI study (did not consent, could not be examined with 
MRI or did not participate for other reasons). An additional 20 patients 
died before acquisition of follow-up MRI scans. Longitudinal MRI data 
was thus available for 87 TBI patients. Only patients with at least 2 suc- 
cessful segmentations of the Tl -weighted 3D brain volumes were in- 
cluded in the study, therefore another 25 patients were excluded. 

For each Tl -weighted 3D scan the segmentation quality was 
assessed visually and scans were excluded if hematomas were included 
as brain tissue, or GM, WM or ventricles were incorrectly delineated. 

See Fig. 1 for details about patient inclusion and exclusion and Fig. 2 
for an example of successful NeuroQuant segmentation. Six scans were 
excluded due to insufficient segmentation quality. Five of these were ac- 
quired in the early phase, and segmentation failed due to hematomas 
segmented as brain tissue (n = 3), incorrect delineation of WM and 
GM of injured frontal lobe (n = 1), and excessive motion (n = 1). 
One 3-month scan failed due to an MRI technical mistake. 

Unsuccessful segmentations resulted in exclusion of two TBI pa- 
tients. Furthermore, in four TBI patients the number of scans was re- 
duced from three to two. 

In total 161 segmented Tl -weighted 3D scans from 62 moderate to 
severe TBI patients were included in the final sample. 37 patients had 
three and 25 patients had two scans successfully segmented. 



2.2. MRI and image processing 

Scanning was performed at three different time points: early phase 
(0-26 days postinjury) and 3 months and 12 months postinjury. The 
early phase scan was acquired as soon as the clinical condition of the pa- 
tients allowed it. 



124 TBI patients 14-65 years old with GCS<13 
Included between October 10, 2004 and October 10, 2007 




20 died (19 died from the head injury, 1 of other 
reason) 

17 did not participate in MRI study (did not consent, 
could not be examined by MRI or did not participate 
for other reasons) 



87 patients participated 



25 excluded 

• 22 because no or only one T 1 W 
NeuroQuant compatible sequence was 
available 

• 1 because of processing failure in 
NeuroQuant 

• 1 because of low quality of the scan 
segmentation, and NeuroQuant 
incompatibility of the other scan 

• 1 due to low quality scan segmentation of 
one scan, and no scan acquired at other time 
point 



62 patients ( 1 7 women) 

• 25 having 2 good quality segmentations included 

• 37 having 3 good quality segmentations included 




Fig. 1. TBI patients with GCS < 13 that could not be explained by other factors than trau- 
matic head injury, admitted to the Neurosurgical Department, St. Olav's Hospital, Trond- 
heim University Hospital, Norway in the period 2004-2007. Flow chart of the inclusion 
and exclusion criteria, and the final number of patients participating in the current study. 



Fig. 2. Successful NeuroQuant segmentation of repeated MRIs obtained in a 21 -year-old 
patient with moderate traumatic brain injury at a) 9 days postinjury and b) 3 months 
and c) 12 months postinjury. 



2.2.1. Scan protocol 

MRI acquisition was performed on a 1.5 T Siemens Symphony Sona- 
ta scanner (Siemens, Erlangen, Germany) with an eight-channel head 
coil. Imaging was performed following a standardized protocol, where 
slice positioning for transversal scans was parallel with the anterior- 
posterior commissural line. A Tl -weighted 3D magnetization prepared 
rapid acquisition gradient echo (MPRAGE) sequence was obtained in 
the sagittal plane with TR = 7.1 ms; TE = 3.45 ms; flip angle = 7°; 
Tl = 1000 ms; FOV = 256 x 256; acquisition matrix of 
256 x 192 x 128, reconstructed to 256 x 256 x 128, giving a recon- 
structed voxel resolution of 1.00 x 1.00 x 1.33 mm. For diffusion- 
weighted imaging a single-shot spin echo planar imaging sequence 
with 19 slices, slice thickness 5 mm, TR = 3300 ms; TE = 110 ms; 
NEX = 4; acquisition matrix size 256 x 256; FoV = 230; in-plane reso- 
lution 0.9 x 0.9 mm and baseline images (b = 0 s/mm 2 ) plus varying 
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diffusion gradient strength along each of three orthogonal directions 
with b = 500 and 1000 s/mm 2 was used. Diffusion trace maps were 
computed from the isotropic diffusion images and used to estimate 
ADC values. FLAIR, T2, and T2*-weighted hemosequences were also 
acquired. 

2.2.2. Segmentation of the Tl -weighted 3D brain scans 

Fully automated segmentation was performed using NeuroQuant 
(CorTechs Labs, La Jolla, CA, US) which is a FDA-cleared tool for clinical 
evaluation of hippocampal atrophy in mild cognitive impairment and 
Alzheimer's disease (Brewer et al., 2012). The output from NeuroQuant 
segmentation contains volumes of the total lobar WM, total cortical GM, 
lateral, third, fourth, and inferior lateral ventricles, cerebellum, hippo- 
campus, amygdala, caudate, putamen, pallidum, thalamus and 
brainstem given as both absolute volume in cm 3 and relative to the in- 
tracranial volume (ICV). In this study only volumes relative to ICV 
(%ICV) are used. The ICV-corrected volumes of the segmented struc- 
tures from the left and right hemispheres were combined. The ICVs in 
the early phase and 12-month scan were significantly different (linear 
mixed model; mean ± standard error: early phase 1740.06 ± 
19.98 cm 3 ; at 12 months 1706.10 ± 19.74 cm 3 ; p < 0.001). ICV is esti- 
mated from the total brain volume in NeuroQuant, and may therefore 
be influenced by brain edema in the early phase and brain atrophy in 
the chronic phase. Differences between ICV measurements based on 
total brain volume in the early and chronic phases are not surprising 
in TBI patients (Ding et al, 2008; Sidaros et al, 2009). Therefore, vol- 
umes of segmented structures were normalized and expressed as pro- 
portions of the ICV estimated from the 3 -month scan. 

Because of the large number of segmented brain structures, we 
chose to analyze the brain structures with low frequency of segmenta- 
tion mistakes. The ICV-corrected volumes of cortical GM, hippocampus, 
brainstem, lenticular nucleus, lobar WM, and ventricles were compared 
between the early, 3-month and 12-month scans. Putamen and 
pallidum were evaluated together as they form the lenticular nucleus. 
The lenticular nucleus was used since it was better delineated than 
the thalamus in the TBI patients. Total ventricular volume was calculat- 
ed as the sum of third, fourth and lateral ventricular volumes. We eval- 
uated the ventricular volumes together since all ventricles have been 
shown to dilate in TBI (Warner et al., 2010b; Kim et al, 2008). 

In addition the volume change over time was calculated for all struc- 
tures. For the solid structures the following was used: early phase vol- 
ume minus volume at 3 months, and early phase volume minus 
volume at 12 months. For the ventricles: volume at 3 months minus 
early phase volume, and volume at 12 months minus early phase vol- 
ume. This was done in order for all mean changes to be positive. 

Total brain atrophy from the early phase to 3 and 1 2 months was es- 
timated using linear mixed model based on the sum of the volume 
changes of the solid structures (cortical GM, hippocampus, lenticular 
nucleus, lobar WM, and brainstem) subtracting the ventricular volume 
change (as ventricular volume change is defined as reversed, please 
see above). 

2.2.3. ADC measurements in GM 

ADC measurements were performed in PACS using Spectra Worksta- 
tion IDS5 1 1.4.1 in all 62 patients. Ten circle-shaped ROIs with radius of 
2.7 mm were positioned in apparently healthy GM bilaterally in the cor- 
tex, hippocampus and the brainstem (Fig. 3). These ROIs were chosen 
because they are located to regions where ADC value changes have 
been reported earlier (see Introduction), the regions are readily identi- 
fiable, and the ROIs are located to parts of GM with minimal WM con- 
tamination. All ROIs were positioned in apparently healthy tissue 
within the structure of interest. In patients with visible focal pathology 
on MRI the ROI was moved to the closest healthy tissue if necessary. 
The cortical ROIs in the superior frontal sulcus and postcentral sulcus 
were positioned at the sulcal depth in order to avoid measuring ADC 
in GM close to the skull with possible minor contusions. The cortical 




Fig. 3. ROI placement for ADC measurements in a) Q superior frontal sulcus ROI, © 
postcentral sulcus ROI; b) © insular gyrus ROI; c) © hippocampal ROI, 0 brainstem ROI. 
Circle-shaped ROIs with radius of 2.7 mm were positioned in apparently healthy tissue. 
In TBI patients with visible focal pathology, the ROI was moved to the closest apparently 
healthy GM tissue within the structure of interest. 

ROIs were placed on the same slice, just above the most dorsal slice 
displaying the lateral ventricles. In the posterior insular gyrus ROIs 
were set on the slice where the head of the caudate nucleus had its 
greatest extent. ROIs in the hippocampus and brainstem were posi- 
tioned on the most caudal slice where mesencephalon was still visible. 
In the hippocampus the ROIs were located as laterally as possible in 
the anterior part of the hippocampal body in order to avoid placing 
the ROIs within the hippocampal head which frequently suffers from 
contusions in TBI. For all ROIs, the ADC values from the left and right 
hemispheres were averaged. The image analyst (V.B.) was blinded to 
clinical information. Intra- rater reliability with regard to ROI positioning 
and subsequent ADC measurements was assessed in 14 scans which 
were reanalyzed after >1 month from the first assessment. Intraclass 
correlation coefficients (ICC) using two-way random single measures 
were used to calculate reliability. 

2.3. Injury -related variables 

GCS score was recorded at hospital admittance, or before intubation 
in case of a pre-hospital intubation. All patients included had GCS < 13 
that could not be explained by other factors than head injury, as ex- 
plained in previous publications from this cohort (Skandsen et al, 
2010; Moen et al., 2012). The head injury severity scale (HISS) is 
based on GCS, where GCS < 8 is considered severe while GCS 9-13 as 
moderate head injury (Stein and Spettell, 1995). 

GCS scores were available for all included patients. 

PTA was recorded and dichotomized into PTA > 2 weeks and 
PTA < 2 weeks which corresponds to moderate and moderate-severe 
Mississippi PTA intervals, respectively. Mississippi PTA intervals have 
been shown to be more accurate predictors of TBI outcome than the tra- 
ditional Russell intervals (Nakase-Richardson et al., 2011). Data on PTA 
were missing in 4 patients. These were excluded from analyses includ- 
ing PTA. 

DAI classification was performed by experienced senior neuroradiol- 
ogists based on the early FLAIR, diffusion-, T2- and T2*-weighted scans, 
for details see Skandsen et al. (2010). DAI was classified into grade 1; 
traumatic lesions confined to lobar WM, grade 2; lesions also detected 
in the corpus callosum, and in grade 3; presence of brainstem lesions 
(Gentry, 1994). DAI classification was available for all included patients. 

2.4. Outcome assessment 

Global outcome was assessed by telephone interview 12 months 
after injury using the structured interview for Glasgow Outcome Scale 
-Extended (GOSE) (Jennettand Bond, 1975; Jennett etal, 1981). To re- 
duce the potential error associated with the telephone setting (LeGrand 
et al., 2007), relatives or caregivers also provided information, in two- 
thirds of the cases, and the best source of information was used 
(based on the judgment of the authors) for details see Skandsen et al. 
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(2010). Outcome data were available for all included individuals, but 
missing in 5 non-included patients. Patients were dichotomized into 
disability (GOSE < 7) and good recovery (GOSE > 7) after 12 months. 

2.5. Statistical analyses 

2.5.1. Characteristics of included and excluded TBI patients. 

The Mann-Whitney If-test, x 2 -test, independent samples T-test and 
the Wilcoxon If-test were used to compare demographic and clinical 
data. 

2.5.2. Longitudinal brain volume changes 

Significant changes in the volume of each of the brain structures at 
each time point was calculated with a linear mixed effects model. 

The relationships between ventricular volume change and volume 
changes of lobar WM, cortical GM and brainstem were explored using 
linear regression separately for the periods early to 3 months and 
early to 12 months. There were separate regression models for lobar 
WM, cortical GM, and brainstem volume changes as the volume chang- 
es were expected to be correlated. 

Impact of TBI subgroup on brain volumes was investigated in the fol- 
lowing dichotomized groups; moderate vs. severe TBI, PTA < 2 weeks 
vs. > 2 weeks, and DAI vs. non-DAI. The volumes were compared over 
time within each group and between the dichotomized groups. For 
DAI a subdivision into four DAI groups, i.e. non-DAI = 0, and DAI grades 
1 , 2, and 3 were performed and the effect of degree of DAI was assessed 
with one-way ANOVA. 

A significant F- value was followed by a paired Student's T-tests with- 
in the group, and with independent samples T-test between the dichot- 
omized groups. The resulting p-values were Bonferroni corrected for 
number of time points evaluated for each structure separately. 

The following analyses of volume differences between the different 
injury mechanism groups were performed: 

2.5.2. 1. Moderate vs. severe TBI. Differences between cortical GM, hippo- 
campal, lenticular nucleus, lobar WM, brainstem and ventricular vol- 
umes in the early phase and at 3 and 12 months were investigated. 
Since DAI was more frequent in severe TBI, the presence of DAI was in- 
cluded in the analysis of lobar WM volume differences between moder- 
ate and severe TBI. 

2.5.2.2. Moderate vs. severe TBI. Differences between cortical GM, hippo- 
campal, and ventricular volumes in the early phase and at 3 and 
12 months were investigated. 

2.5.2.3. DAI vs. non-DAI. Differences in cortical GM, lobar WM, brainstem, 
and ventricular volumes in the early phase and at 3 and 1 2 months were 
investigated. 

2.5.2.4. DAI grades 0-3. The effect of degree of DAI (no DAI, DAI grade 1, 2 
or 3) as an ordinal variable on brainstem and total ventricular volumes 
was assessed in the early phase and at 3 and 12 months. 

2.5.3. Differences in volume changes over time between brain structures 
The presence of statistically significant differences in rate of volume 

loss between different brain structures was evaluated in a linear mixed 
model with random intercept since this model fitted the data best, i.e. 
had the lowest Akaike's information criterion. There were not enough 
observations to build a model exploring non-linear changes over time. 
The following comparisons were done: cortical GM vs. hippocampus, 
cortical GM vs. lobar WM, hippocampus vs. lentiform nucleus, and hip- 
pocampus vs. lobar WM matter. For these analyses significant differ- 
ences in the interaction term between volume loss over time and 
brain structure were calculated, and the estimated difference in slope 
with 95% confidence intervals was reported. 



2.5.4. Outcome prediction 

Hierarchical regression analysis was performed to determine if 3- 
month MRI-derived volumes contributed significantly to 12-month 
GOSE outcome after taking into account the established outcome pre- 
dictors GCS score, PTA (PTA < 2 weeks = 0, PTA > 2 weeks = 1 ) and 
age which were entered in the first step as the base model with which 
all subsequent models were compared. The MRI measures were cortical 
GM, hippocampal, lenticular nucleus and ventricular volumes at 
3 months. In addition we used the Rotterdam CT scores (Maas et al, 
2005) obtained from the worst CT scan, to see if this neuroimaging mo- 
dality provided additional predictive value (Yuh et al, 2012). Overall 
variance explained {R 2 ) and statistical significance were calculated for 
every model. K 2 -change and significance of F-value change were calcu- 
lated as differences between the base model and all subsequent models. 

2.5.5. Longitudinal analysis of ADC changes in GM 

Intraclass correlation coefficients (ICC) using two way random sin- 
gle measures were used to calculate reliability of the ROI ADC values 
(Gwet, 2008). 

A linear mixed effects model was used in the statistical comparisons 
of the ADC values in the different ROIs over time in the entire TBI group, 
and between the TBI subgroups; moderate vs. severe TBI, PTA < 2 weeks 
vs. >2 weeks, and DAI vs. non-DAI. The effect of degree of DAI (non- 
DAI = 0, DAI grades 1, 2, and 3) was assessed with ANOVA. Significant 
F-values were followed by paired Student's T-test or independent sam- 
ple T-test. The resulting p-values were Bonferroni corrected as described 
earlier. Associations between ADC values and corresponding volumes at 
3 and 12 months were examined using Pearson correlation. 

A similar hierarchical regression analysis as described for the vol- 
umes above was performed to establish if cortical ADC values from the 
superior frontal and post-central sulci ROIs at 3 months were indepen- 
dent significant predictors of 12 month outcome. 

Clinical and demographic data (age, GCS and GOSE scores) are given 
as mean ± standard deviation, and volumes and ADC values as mean ± 
standard error. The threshold for statistical significance was p < 0.05 for 
all analyses, including analyses applying Bonferroni correction for mul- 
tiple comparisons. 

3. Results 

3.1. TBI patient characteristics (Table 1 ) 

There were no significant differences with regard to the clinical var- 
iables between included and excluded TBI patients who survived for 
12 months. Of the included patients 58% suffered moderate and 42% se- 
vere TBI. See Table 1 for details. 

Comparison between the included males (n = 45) and females (n = 
17) showed that GCS score, HISS, and PTA and GOSE score at 12 months 
were similar between the sexes. 

Injury mechanism, age and number of days between the head injury 
and the MRI scans were similar in the moderate and severe TBI patients. 
There were significantly more patients with PTA > 2 weeks, any grade of 
DAI and poor 1 2-month outcome in the severe than in the moderate TBI 
group. 

3.2. Longitudinal brain structure volume changes in the entire TBI group 
(Figs. 4 and 5) 

During the first 12 months after TBI significant volume differences 
were found for all analyzed structures (Fig. 4). Cortical GM experienced 
a rapid volume reduction and was significantly reduced at 3 months 
compared to the early scan, but did not change significantly after that. 
The hippocampal and lenticular nucleus volumes decreased more slow- 
ly and the volume was not significantly reduced before 12 months 
postinjury. Lobar WM volume was significantly reduced first at 
12 months compared to both the early and the 3-month scan. Brainstem 
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Table 1 

TBI patient characteristics. 







All (n = 62) 


Moderate a (n = 36) 


Severe' 3 (n — 26) 


Age (years) 




32.5 ± 15.2 


34.3 ± 15.8 


30.1 ± 14.4 


GCS 




9.1 ± 3.4 


11.5 ± 1.6* 


5.7 ± 1.8* 


Mechanism of injury 


Road accident 


34 (55%) 


19(53%) 


15 (58%) 




Fall injury 


23 (37%) 


13 (36%) 


10(38%) 




Other 


4 (6%) 


3 (8%) 


1 (4%) 




Unknown 


1 (2%) 


1(3%) 


0 


Post-traumatic amnesia 


>2 weeks 


19(31%) 


5(14%)* 


14 (54%)* 




<2 weeks 


39 (63%) 


30 (83%)* 


9 (35%)* 


Diffuse axonal injury 


No 


19(31%) 


17(47%)* 


2 (8%)* 




Yes 


41 (66%) 


19(53%)* 


22 (85%)* 




Grade 1 


13 (21%) 


5 (14%)* 


8 (31%)* 




Grade 2 


17(27%) 


11 (31%)* 


6 (23%)* 


Outcome at 12 months 


Grade 3 


11 (18%) 


3 (8%)* 


8 (31%)* 


With disability (GOSE < 7) 




27 (44%) 


12(33%)* 


15 (58%)* 


Good recovery (GOSE 7-8) 




35 (56%) 


24 (67%)* 


11 (42%)* 


Days between injury and MRI C 


Scan 1 


10.2 ± 6.7 


9.0 ± 6.2 


12.1 ± 7.3 




Scan 2 


98.4 ± 15.2 


98.2 ± 14.3 


98.6 ± 16.7 




Scan 3 


373.5 ± 26.3 


377.1 ± 30.1 


369.0 ± 20.3 



a GCS 9-13. 
b GCS 3-8. 
c In days. 

Significant difference between moderate and severe TBI. Numbers of patients (percentages) or means ± standard deviations are shown. Percentages were calculated considering the 
total number of patients in each column. GCS, Glasgow Coma Scale; GOSE, Glasgow Outcome Scale— Extended; MRI, magnetic resonance imaging. 



volumes decreased significantly between the early and 3-month scans 
and also between the 3-month and 12-month scans. The ventricular 
volume expanded from the early phase, but a significant increase was 
present first at 12-months. Statistical comparison of the interaction be- 
tween slope of volume loss over time for the different brain structures 
demonstrated that cortical GM loss followed a significantly different tra- 
jectory compared to the hippocampus (estimated difference in slope 
was 1.2 (CI: 0.7-1.8), p < 0.0001), but not compared to lobar WM (esti- 
mated difference in slope was —0.04 (CI: —0.9-0.8); p = 0.9). The hip- 
pocampus and lenticular nucleus had similar trajectories of volume loss 
(estimated difference in slope was —0.01 (CI: —0.3-0.0); p = 0.1). For 
lobar WM and hippocampus significantly different trajectories of vol- 
ume loss were present (estimated difference in slope was 1.2 (CI: 
0.7-1.7); p< 0.0001). 

Ventricular dilation was strongly associated with change in 
brainstem volumes at 3 months (R 2 = 0.246, p = 0.001 ) and 12 months 
(R 2 = 0.488, p < 0.001 ), and total cortical GM volume at 3 months (R 2 = 
0.143,p = 0.012) and 12 months {R 2 = 0.331, p< 0.001). Ventricular di- 
lation was not significantly associated with lobar WM volume changes 
at any time point. 

The total brain volume decline was 3.8 ± 0.8% of ICV between the 
early and 3-month scans, and 5.7 ± 0.8% of ICV between the early and 
12-month scans. 



3.3. The impact of clinical TBI subtypes on brain volumes (Figs. 6 and 7) 
3.3.1. TBI severity 

Lobar WM and brainstem volume at 12 months were smaller in the 
severe compared to the moderate TBI group (lobar WM: 26.34 ± 0.58% 
ICV vs. 28.60 ± 0.51% ICV, p = 0.012; brainstem: 1.36 ± 0.03% ICV vs. 
1.47 ± 0.02% ICV, p = 0.005) while cortical GM, ventricular, hippocam- 
pal and lenticular nucleus volumes were similar. The difference in lobar 
WM volume at 12 months was significant also after correcting for the 
presence of DAI (p = 0.012). 

Moreover, hippocampal volume change between the early and 3- 
month scans was larger in the severe group compared to the moderate 
group (0.03 ± 0.01 %ICVvs. 0.00 ± 0.01%ICV,p = 0.046). For the other 
structures, there were no differences in volume changes over time be- 
tween moderate and severe TBI. 



3.3.2. PTA duration (Fig. 6a, b) 

Patients with PTA > 2 weeks had smaller hippocampi at 3 and 
12 months, and larger ventricles at 3 and 12 months compared to pa- 
tients with PTA < 2 weeks. Cortical GM volumes did not differ signifi- 
cantly between the PTA groups. 

Hippocampal volume change over time was larger in the 
PTA > 2 weeks than in the PTA < 2 weeks group between the early 
and 3-month scans (0.05 ± 0.01 %ICV vs. 0.00 ± 0.01 %ICV, p = 
0.003). Cortical GM and ventricular volume changes over time did not 
differ between the PTA groups. 



3.3.3. DAI groups (Figs. 6c, d and 7) 

Smaller brainstem volumes were present at all time points in the DAI 
group compared to the non-DAI group. No differences in cortical GM, 
lobar WM or ventricular volumes were present between the DAI and 
non-DAI groups. 

The non-DAI group had larger brainstem volumes than the DAI 
grade 2 group in the early phase; than the DAI grade 2 and 3 groups at 
3 months, and the DAI grade 1, 2 and 3 groups at 12 months. In DAI 3 
brainstem volumes decreased significantly between the early and 
3 month scans, and further till 12 months. In DAI 2 brainstem volume 
decrease was present first after 12 months. In DAI 1 brainstem volumes 
decreased between the early scan and 12 months and between the 3 
and 12 month scans. In the non-DAI group, brainstem volumes did not 
change with time. 

No significant differences between the DAI groups were observed in 
ventricular volumes. 

For the non-DAI, DAI 2 and DAI 3 groups the ventricular volume de- 
creased between the early and 3 month scans, and between the early 
and 12 month scans. DAI 1 ventricular volume decreased between the 
early and 12 month scans. 

3.4. Outcome prediction based on 3-month MRI volumes (Table 2) 

The hierarchical regression analyses with PTA, GCS scores and age in 
the base model demonstrated that the MRI-derived brain volumes of 
cortical GM, hippocampus, and lenticular nucleus at 3 months explained 
significantly more of the final outcome than the base model (Table 2). 
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Fig. 4. Longitudinal changes in mean ICV-corrected brain volumes with standard errors in 
the early phase (0-26 days post-injury, scan 1) and 3 months (scan 2), and 12 months 
(scan 3) postinjury. Volumes of the left and right sides of the structures were summed 
and evaluated together. ***p < 0.001, **p < 0.01, *p < 0.05. 



The ventricular 3 -month volume did not add to the 12-month outcome 
prediction, and neither did the Rotterdam CT score. 

3.5. Analysis of ADC value changes in GM 

3.5 J. Intra-rater reliability of ADC measurements in the ROIs 

ICC of repeated measurements were 0.74 for superior frontal sulcus, 
0.86 for postcentral sulcus, 0.69 for insular gyrus, 0.88 for hippocampus 
and 0.89 for brainstem. 



3.5.2. Overall changes in entire TBI group 

No ADC value changes over time were found in any ROI in the entire 
TBI group. 

3.5.3. The impact of clinical TBI subtypes on ADC value changes 

3.5.3.1. TBI severity (Table 3a). The severe TBI group had higher ADC 
values in the brainstem at 3 months, the superior frontal sulcus at 3 
and 12 months, and the postcentral sulcus at 12 months compared to 
the moderate TBI group. No significant differences in ADC values be- 
tween moderate and severe TBI groups were found at any time point 
for the insula or hippocampus. 

In the severe TBI group the ADC values increased over time in the su- 
perior frontal sulcus from the early phase to 3 months (p = 0.043), and 
the postcentral sulcus from the early phase to both 3 and 12 months 
(early phase < 3 months, p = 0.028; early phase < 12 months, p = 
0.027). There were no differences in brainstem, insula or hippocampal 
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Fig. 5. Illustration of the volume changes over time. Volume changes (as %ICV) early 
phase-3 months and 3 months-12 months postinjury. The differences between different 
time points are calculated using linear mixed model. Means ± error bars are plotted. 



ADC values between any time points. There were no differences in 
ADC values over time in any ROI in the moderate TBI group. 

3.532. PTA duration (Table 3b). At 12 months the ADC values in the 
superior frontal sulcus were higher in patients with PTA > 2 weeks 
than in those with PTA < 2 weeks, while no differences were found 
in the other ROIs. In the PTA > 2 weeks group ADC values in insula in- 
creased between 3 and 12 months (p = 0.046). In the group with 
PTA < 2 weeks no differences in ADC values overtime were found 
in any ROI. 

3.5.4. DAI groups (Figs. 8 and 9) 

In the superior frontal sulcus and postcentral sulcus ADC values 
were higher in TBI patients with DAI than in non-DAI patients at 3 and 
12 months. The ADC values increased significantly from the early 
phase to 3 months in the superior frontal (p = 0.009) and the 
postcentral sulcus (p = 0.005) in the DAI group. No differences were 
found in ADC values or ADC value change over time in the insula, hippo- 
campus or brainstem ROIs between the non-DAI and the DAI groups, or 
within the DAI grades. 

The different DAI grades were associated with significant differences 
in cortical GM ADC values (Fig. 9). In DAI 3, the superior frontal sulcus 
ADC values were higher than those in the non-DAI group at 3 months, 
and compared to both non-DAI and DAI 1 at 12 months. In DAI 2 the 
ADC values were higher in the superior frontal sulcus at 3 and 1 2 months 
than those in the non-DAI group. The DAI 3 group had higher ADC 
values in the postcentral sulcus at both 3 and 12 months than non- 
DAI and DAI 1 groups. In the superior frontal sulcus ADC values de- 
creased between 3 and 12 months (p = 0.018) in DAI 1, and increased 
between the early phase and both 3 and 12 months in DAI 2 (early 
phase < 3 months, p = 0.036; early phase < 12 months, p = 0.018). 
There were no significant changes over time within the non-DAI and 
the DAI 3 groups. 
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Fig. 6. Mean volumes of different brain structures plotted for: a) hippocampus and b) ven- 
tricles of patients with PTA > and <2 weeks, c) Brainstem and d) ventricles of DAI and non- 
DAI patients. For statistical evaluation of volume differences see Materials and methods 
and Results. 

The ADC values increased between the early phase and 3 months in 
the DAI 2 group (p = 0.020) in the postcentral sulcus. In DAI 3, the ADC 
values were significantly increased both at 3 and 12 months compared 
to the early phase (early phase < 3 months, p = 0.003; early 
phase < 12 months, p = 0.003). ADC values in postcentral sulcus did 
not change with time in the non-DAI and DAI 1 groups. 



3.5.5. ROIADC values vs. brain structure volumes. 

There were no correlations between the ADC values in the cortex, 
hippocampus or brainstem and the volumes of these structures at 3 
and 12 months. 



3.5.6. ADC and outcome prediction 

ADC values of superior frontal and postcentral sulci from 3 months 
did not add to the 12-month outcome prediction (Table 4). 
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Fig. 7. The association between the longitudinal changes in brainstem and ventricular ICV- 
corrected volumes for the non-DAI and the three DAI groups (DAI grades 1, 2, and 3). 
Figures show means and standard errors. Left and right side volumes were summed and 
evaluated together. Significant differences in brainstem and ventricular volumes between 
the non-DAI and the three DAI groups were found at all time points. Please note that the 
units on the y axis differ between the plots, and that the x axis does not cross the y axis 
at 0. 



Table 2 

Hierarchical regression analysis of 12-month outcome prediction using the 3-month brain 
structure volumes with GCS and PTA as base model. 



Models 


R 2 


K 2 -change 


Significance of F change 


Base model 


0.385 






CortGM 3 months 


0.515 


0.053 


0.023 


Ventricles 3 months 


0.467 


0.005 


0.494 


Hippocampus 3 months 


0.535 


0.073 


0.007 


Lenticular ncl. 3 months 


0.529 


0.067 


0.010 


Rotterdam CT score 


0.465 


0.003 


0.596 



Base predictors: GCS, PTA duration (0 = PTA < 2 weeks; 1 = PTA > 2 weeks), and age 
CortGM predictors: GCS, PTA duration, age, and cortical GM volume at 3 months 
Ventricular predictors: GCS, PTA duration, age, and ventricular volume at 3 months 
Hippocampal predictors: GCS, PTA duration, age, and hippocampal volume at 3 months 
Lenticular nucleus predictors: GCS, PTA duration, age, and lenticular nucleus volume at 
3 months 

Rotterdam CT score predictors: GCS, PTA duration, age, and Rotterdam CT score 
CortGM, cortical gray matter; GCS, Glasgow Coma Scale; PTA, posttraumatic amnesia du- 
ration. 



4. Discussion 

This is to the best of our knowledge the largest number of moderate 
and severe TBI patients included in a prospective, longitudinal MRI 
study of brain volumes and diffusion changes where data have been sys- 
tematically collected at 3 different time points, i.e. early phase and 3 and 
12 months postinjury. Moreover, the patient group included in this 
study was representative for the entire group of TBI patients surviving 
for 1 year. The results not only confirmed most of our predictions, but 
also revealed new, unexpected findings. 



4.1. Longitudinal volume changes 

Significant cortical GM volume loss took place between the early and 
3-month scans without any further significant decline at 12 months. 
This is the first direct evidence for cortical GM volume loss being an 
early occurrence after TBI in humans, similar to results in animal models 
(Rodriguez-Paez et al, 2005; Smith et al, 1997; Liu et al, 2010). The 
rapid cortical volume loss probably stems from a combination of neuro- 
nal necrosis and apoptosis initiated at time of injury and peaking early 
postinjury (Conti et al, 1998). Supporting these observations in ani- 
mals, human PET studies show early and widespread neurodegenera- 
tion in intact cortex following TBI ( Wu et al, 2004; Bergsneider et al., 
2001 ). We did not find any evidence for GCS scores, duration of PTA 
or presence of DAI being associated with cortical volume loss. The effect 
of DAI and PTA on cortical volumes may, however, be localized to specif- 
ic cortical regions for instance connected to affected WM tracts or the 
medial temporal lobe, and hence not be reflected in total cortical 
volumes. 

The volume loss in the hippocampus was significantly slower than 
that in the cortex and not significantly different from the early and 
3 month scans before 12 months postinjury. This slow, continuous 
loss of hippocampal volume after TBI in humans concurs with animal 
studies (Immonen et al., 2009; Smith et al, 1997). However, both ani- 
mal data (Immonen et al, 2009; Smith et al., 1997) and autopsy data 
from humans (Maxwell et al., 2003) demonstrate significant early neu- 
ronal loss in the hippocampus. The lack of a significant decline in hippo- 
campal volume between the early and 3 month-scans in the current 
study may be due to the early scans being obtained at -10.2 days 
after injury. At that time point some volume loss may already have 
taken place. Alternatively, MRI-based hippocampal volume loss may 
take longer to develop than histological identified cell loss. Still, our re- 
sult concurs with findings in cross-sectional studies where decreased 
hippocampal volume is detected first after 3 months in moderate-se- 
vere TBI patients compared to matched controls (Wilde et al, 2007; 
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Table 3 

Longitudinal ADC values of TBI subgroups and comparison between the clinical TBI subtypes. 



a) Moderate vs. severe TBI. 







Moderate TBI 


Severe TBI 


p-Value 


Superior rrontal sulcus 




















Early phase 


68.2 


± 


1.1 


68.1 


± 


1.3 


MC 

IN J 




mnnthc 


66.9 


± 


1.0 


71.5 


± 


1.2 


0.011 




1 2 months 


65.0 


± 


1.1 


70.7 


± 


1.2 


0 005 


Postcentral sulcus 




















Early phase 


70.2 


± 


0.9 


69.6 


± 


1.1 


NS 




3 months 


69.9 


± 


0.9 


73.2 


± 


1.0 


IN j 




1 2 months 


68.9 


± 


1.0 


73.2 


± 


1.0 


0 009 


Insular gyrus 




















Early phase 


81.5 


± 


1.1 


81.9 


± 


1.3 


NS 




3 months 


80.7 


± 


1.0 


81.9 


± 


1.2 


NS 




12 months 


81.4 


± 


1.1 


83.9 


± 


1.2 


NS 


Hippocampus 




















Early phase 


86.1 


± 


1.4 


88.0 


± 


1.7 


NS 




3 months 


85.8 


± 


1.3 


87.8 


± 


1.5 


NS 




12 months 


85.1 


± 


1.4 


88.7 


± 


1.5 


NS 


Brainstem 




















Early phase 


80.6 


± 


1.7 


81.5 


± 


2.1 


NS 




3 months 


80.5 


± 


1.6 


87.5 


± 


1.9 


0.016 




12 months 


80.3 


± 


1.7 


83.5 


± 


1.9 


NS 


b) Patients with PTA < 2 weeks vs. >2 weeks 






PTA < 2 weeks 


PTA > 2 weeks 


p-Value 


Superior frontal sulcus 




















Early phase 


68.2 


± 


1.0 


67.9 


± 


1.5 


NS 




3 months 


68.0 


± 


1.0 


70.4 


± 


1.5 


NS 




12 months 


66.3 


± 


1.1 


70.8 


± 


1.5 


0.037 


Postcentral sulcus 




















Early phase 


70.0 


± 


0.9 


69.8 


± 


1.4 


NS 




3 months 


70.4 


± 


0.9 


72.9 


± 


1.3 


NS 




12 months 


69.9 


± 


1.0 


72.6 


± 


1.3 


NS 


Insular gyrus 




















Early phase 


81.7 


± 


0.9 


81.5 


± 


1.4 


NS 




3 months 


81.9 


± 


0.9 


79.5 


± 


1.3 


NS 




12 months 


81.8 


± 


1.0 


84.0 


± 


1.3 


NS 


Hippocampus 




















Early phase 


87.0 


± 


1.1 


86.2 


± 


1.6 


NS 




3 months 


86.0 


± 


1.1 


85.0 


± 


1.5 


NS 




12 months 


86.1 


± 


1.1 


87.8 


± 


1.5 


NS 


Brainstem 




















Early phase 


80.8 


± 


1.5 


81.3 


± 


2.2 


NS 




3 months 


81.8 


± 


1.4 


83.8 


± 


2.1 


NS 




12 months 


81.2 


± 


1.5 


82.6 


± 


2.1 


NS 



ADC in 10~ 5 m 2 /s. Values given as means ± standard deviations. 
Early phase scans were acquired 0-26 days postinjury. 

Statistical differences between moderate vs. severe or PTA < 2 vs. >2 weeks were Bonferroni-corrected for every structure separately. 
ADC, apparent diffusion coefficient; NS, non-significant; TBI, traumatic brain injury. 



Bigler et al, 1997). The current data also agree with results from longi- 
tudinal studies with two time points, which show that hippocampal vol- 
umes decline over a protracted period from 1 week to 2.5 years (Warner 
et al., 2010b; Ng et al., 2008). The slowly evolving hippocampal volume 
loss is consistent with the protracted apoptotic neuronal death de- 
scribed in this region (Williams et al., 2001 ), coupled with glial prolifer- 
ation (Gao and Chen, 2013) which may counteract some of the volume 
change, and impaired neurogenesis (Gao et al., 2008) which may affect 
volume when accumulated over time. The hippocampus is critical for 
declarative memory (Squire et al, 2004; Eichenbaum et al, 2007), and 
hippocampal damage is known to produce symptoms of anterograde 
and retrograde amnesia (Cipolotti and Bird, 2006). It was thus not sur- 
prising to find that patients with PTA > 2 weeks had significantly 
lower hippocampal volumes at both 3 and 12 months than patients 
with PTA < 2 weeks. Moreover, hippocampal volume loss was markedly 
more pronounced from the early to 3 month scans in the PTA > 2 weeks 
group. These results demonstrated that PTA reflects hippocampal neu- 
ronal injury and that more severe PTA is associated with more rapid 
and extensive neuronal death in this region. 



Significant volume loss developed over 12 months in the lenticular 
nucleus. Decreased lenticular nucleus volume has been reported in 
cross-sectional and longitudinal studies of TBI in children and adults 
>1 year after injury (Wilde et al, 2007; Anderson and Bigler, 1995; 
Primus et al., 1997). The histopathological changes in the lenticular nu- 
cleus include apoptosis, neuronal loss and/or shrinkage due to 
transneuronal degeneration and loss of myelinated axons as described 
in animal studies (Tong et al, 2002; Ng et al, 1994). Based on the cur- 
rent results these processes take place over a protracted period. The 
present data do not support the notion that basal ganglia respond 
more rapidly and briefly to TBI than the hippocampus (Hicks et al., 
1996). Rather, the current study showed that the trajectory of volume 
loss was similar in the subcortical GM structures, i.e. the lenticular nu- 
cleus and hippocampus. Moreover, the subcortical GM structures were 
shown to differ significantly from the cortex with regard to the slopes 
of volume loss suggesting important differences in the histopathological 
responses between these GM structures to TBI. 

As predicted lobar WM volume loss decreased mainly in the late 
phase, being significant first at 12 months. This concurs with the slowly 
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Fig. 8. Longitudinal changes in mean ADC values with standard deviations in non-DAI (□) 
and DAI (■) patients in the superior frontal sulcus (a), postcentral sulcus (b), insular gyrus 
(c), hippocampus (d) and brainstem (e). Right and left side values were averaged. Statis- 
tical significant differences between non-DAI and DAI groups are marked with *. Signifi- 
cant differences were found between: a) DAI (1) and DAI (2), p = 0.009; non-DAI 
(2) and DAI (2),p = 0.028; non-DAI (3) and DAI (3),p = 0.001. b) DAI (1) and DAI (2), 
p = 0.006; non-DAI (2) and DAI (2), p = 0.020; non-DAI (3) and DAI (3), p = 0.004. 
The p-values were corrected for multiple comparisons. Please note that the units on the 
y axis differ between the plots, and that the x axis does not cross they axis at 0. 



evolving pathological response in WM (Vargas and Barres, 2007; 
Coleman, 2005; Coleman and Freeman, 2010). The extended period of 
volume loss found in lobar WM in humans in the current study is 
most likely the macroscopic consequence of the protracted histological 
changes in WM. Previous longitudinal studies have only been able to 
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Fig. 9. The associations between longitudinal changes in mean ADC values for the non-DAI 
group (non-DAI = DAI 0) and the three DAI groups (DAI grades 1, 2, 3). Longitudinal anal- 
ysis within each DAI group was performed separately for both ROIs using a linear mixed 
model. Significant F-values were followed by the Student's T-test. The resulting p-values 
were Bonferroni corrected separately for every ROI. At each time point, differences be- 
tween DAI subgroups were evaluated using one-way ANOVA, and significant F-values 
were followed by independent samples T-tests, the resulting p-values were Bonferroni 
corrected for the three time points separately. Significant differences between different 
time points for each DAI group were confirmed and marked as: ©—significantly different 
from scan 2; $— significantly different from scan 3 (***p < 0.001, **p < 0.01, *p < 0.05). 



Table 4 

Hierarchical regression analysis of 12-month outcome prediction using the 3-month ADC 
values with GCS, PTA and age as base model. 



Models 


R 2 


K 2 -change 


Significance of F change 


Base model 


0.462 






Superior frontal sulcus model 


0.479 


0.017 


0.208 


Postcentral sulcus model 


0.480 


0.018 


0.193 



Base predictors: GCS, PTA duration, and age. 

Superior frontal sulcus model: GCS, PTA duration, age, and superior frontal sulcus ADC 
value at 3 months 

Postcentral sulcus model: GCS, PTA duration, age, and postcentral sulcus ADC value at 
3 months 

ADC, apparent diffusion coefficient; GCS, Glasgow Coma Scale; PTA, posttraumatic amne- 
sia duration. 



show a significant decline in WM volumes between an early (1 — 

2 months) and a later time point (6-12 months) after mild to severe 
TBI (Ding et al, 2008; Sidaros et al, 2009; Warner et al., 2010b). Surpris- 
ingly the trajectory of WM loss and cortical GM loss did not differ signif- 
icantly although the slope of WM loss was significantly different from 
the hippocampus, lenticular nucleus, and brainstem. This finding 
could reflect a significant interaction between WM volume loss and cor- 
tical volume loss over time at the individual level. Another explanation 
could be that volume losses in cortical GM and lobar WM volumes were 
the least linear. The lack of a significant interaction between volume 
losses over time between these two regions could hence be due to the 
model not being able to take nonlinearities into account. 

Lobar WM volume loss was associated with injury severity as pa- 
tients with severe TBI had significantly reduced lobar WM volumes 
compared to the moderate TBI group. This difference was present also 
after correcting for the presence of DAI. Unexpectedly, DAI was not as- 
sociated with decreased lobar WM volume. This may be due to DAI clas- 
sifications being based on visible lesions on conventional MRI which 
may not fully describe WM involvement in TBI. 

The brainstem was the only structure in which volume loss was sig- 
nificant between every scan session, thus demonstrating particularly 
pervasive and enduring consequences of TBI to this region. Previous 
MRI studies have reported brainstem volume loss in cross-sectional 
and longitudinal studies with two time points in both children 
and adults with moderate and/or severe TBI (Sidaros et al., 2009; 
Warner et al, 2010b; Fearing et al., 2008). WM injury probably con- 
tributed to the early brainstem volume decline since patients with 
DAI grade 3 had the fastest and largest reduction of brainstem vol- 
ume, present from 3 months. Also the group with severe TBI had 
smaller brainstem volumes, but this was detectable first after 
12 months, appearing at the same time as reduced lobar WM vol- 
umes in this group. The WM tracts in the brain stem thus appear 
to be particularly sensitive to TBI. 

Ventricular volumes were significantly increased at 3 months, 
followed by a non-significant increase at 12 months. This result is a re- 
finement of the previously described time frame of significant ventricu- 
lar dilation which was shown to last until ~2-7 months postinjury 
(Blatter et al., 1997). Ventricular dilation at 3 and 12 months was asso- 
ciated with brainstem and cortical GM volumes. Neither lobar WM vol- 
umes nor the presence of DAI was associated with larger ventricular 
volumes. These results demonstrate the importance of deep WM loss 
in particular, followed by cortical GM to ventricular dilation in TBI, rath- 
er than hemispheric WM loss. 

For the entire TBI group mean brain atrophy was -4% of ICV after 

3 months and -6% of ICV after 12 months, demonstrating that the bulk 
volume loss occurred within the first 3 months. Still, a notable 
protracted component to the atrophy was present. No directly compara- 
ble human studies exist, but endpoint brain volume atrophy in the 
chronic phase after TBI is reported to be from 1.4% to 8% in cross- 
sectional and longitudinal studies in mild to severe TBI (Trivedi et al., 
2007; Ding et al., 2008; Sidaros et al, 2009). 
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4.2. longitudinal ADC changes 

ADC values increased over time in the individuals with the most se- 
rious brain injuries, i.e. severe TBI, PTA > 2 weeks, and DAI, and there 
were no associations between ADC values and volumes. Participants 
with severe TBI had higher ADC values in deep, radiologically unaffected 
sulcal cortical GM, which increased significantly throughout the 
12 month period, even when total cortical volume remained similar be- 
tween 3 and 12 months. Thus it appears that small-scale remodeling of 
cortex after TBI persisted for a longer period than volumetric changes. 
Moreover, a significant contribution from DAI to the increase in cortical 
ADC values was present. More severe DAI grades gave rise to increased 
cortical ADC values developing at a more rapid timescale. In DAI the cor- 
tical ADC changes may be connected to delayed neuronal death, soma 
shrinkage, loss of cortical axons and myelination developing over an ex- 
tended time period in this condition (Goetz et al., 2004; Hergan et al, 
2002). 

As expected there were increased ADC values at 3 months in the 
brainstem in severe TBI, but not in moderate TBI. This finding lends fur- 
ther support to the notion that deep brain injury is significantly more 
pronounced in severe than moderate TBI. PTA > 2 weeks was not asso- 
ciated with increased hippocampal ADC values. Both increased and de- 
creased hippocampal ADC values are reported after TBI in animals and 
humans (Bendlin et al., 2008; Obenaus et al., 2007; Van putten et al., 
2005). We predicted that PTA > 2 weeks would lead to increased ADC 
values due to increased cell loss, but did not find support for this as- 
sumption in the data. It is possible that the pathophysiological mecha- 
nisms leading to PTA are more heterogeneous than we expected and/ 
or that volume loss combined with glial proliferation, inflammation 
(Woodcock and Morganti-Kossmann, 2013; Lin and Wen, 2013; 
Giunta et al., 2012) and/or neurogenesis (Gao et al, 2008; Greer et al., 
2011; Gao and Chen, 2013) counteracted any effect of neuronal loss 
on ADC values. 

4.3. MRI-derived 3-month volumes for 12-month outcome prediction 

The hierarchical regression analysis demonstrated that MRI-derived 
volumes at 3 months contributed independently to 12-month outcome 
prediction even adjusting for duration of PTA, GCS scores and age. This 
finding demonstrates a significant clinical potential for MRI volumes ob- 
tained at 3 months for outcome prediction. Rotterdam CT score, on the 
other hand, did not have additional predictive value. Furthermore, cor- 
tical ADC values at 3 months did not improve outcome prediction. 

4.4. Methodological considerations 

This is the first study using a fully automated method, NeuroQuant, 
for deriving brain volumes from MRI scans in patients with TBI. Our re- 
sults show that NeuroQuant performed very well in moderate-severe 
TBI patients even in the presence of different types of injuries and post- 
operative changes. Its limitations were restricted to segmentation prob- 
lems in the early phase due to hematomas and contusions, and 
excessive movement during scanning. For the 3-month scans no seg- 
mentations failed. Since several volumes derived from the 3-month 
scan were independent predictors of outcome at 12 months, MRI at 
3 months appears to be of clinical value in the assessment of the brain 
injury and prediction of later function. However, for the individual pa- 
tient the presence of focal injury (contusions, hemorrhages) also plays 
a significant role for outcome, but it was beyond the scope of the present 
paper to evaluate the interaction between focal injuries and the atrophy 
in different brain structures. 

It should be noted that the volume change measures obtained using 
independent NeuroQuant segmentations would not be expected to 
have the spatial specificity and power to detect subtle change that 
many across time point registration methods might provide (Holland 
et al, 2009; Fox et al, 2001; Thompson et al, 2004). On the other 



hand, NeuroQuant has previously been shown to have results compara- 
ble to that of hand segmentation of subregional volumes by anatomical 
experts (Brewer et al, 2009), and thus provides an identifiable and 
translatable measure of structure volume. As a FDA-cleared tool used 
at more than 100 clinical sites and applied to over 30,000 individuals 
in routine clinical practice, the tool yields standardized data that are di- 
rectly comparable to results obtained in the clinic. Further, by providing 
independent segmentations of each time point, this approach is inverse- 
consistent and not subject to registration bias that may inflate measures 
of change provided by registration-based techniques (Holland et al., 
2012). Thus, the approach yields a conservative measure of the regional 
brain volume changes associated with TBI. 

Reliability of the ADC measurements was high to very high, but only 
performed by a single rater and interrater variability could not be 
assessed as Ozturk et al. (2008) recommend. There are several ways of 
analyzing ADC values, and we choose manually positioned ROIs to try 
to optimize ROI localization to avoid partial volume effects as slice thick- 
ness was 5 mm. However, more isotropic voxels and/or smaller voxels 
would be preferable to avoid partial volume effects, which potentially 
could have made it easier to detect statistical significant differences be- 
tween groups and over time. For all the GM ROIs, the effect sizes of the 
differences between ADC values for the different TBI groups were very 
large. This was also seen in the regions where there were no statistical 
significant differences. Taken together, our results show that that vari- 
ability in ADC measures is substantial in the material, and that ROI posi- 
tioning probably plays some part in this, but that the largest 
contribution comes from individual differences within and between 
subjects. ADC values in GM may change slightly during normal aging. 
However, according to earlier studies, brain ADC values of healthy sub- 
jects remain constant during the age period 20-60 years (Klimas et al., 
2013; Watanabe et al, 2013). As most of the subjects included in the 
present study are in this age range, we assumed that the diffusion 
changes revealed are TBI-associated. Including a matched control 
group may have helped in the interpretation of the ADC values, which 
relies on inferences from histopathology. Newer diffusion based MRI 
methods such as restriction spectrum imaging (White et al., 2013) 
may be better suited for characterizing tissue microstructural changes 
in TBI in the future. 

The large number of analyses conducted across multiple brain re- 
gions increased the susceptibility to type I errors in the current study. 



5. Conclusion 

In summary, the effect of TBI on brain volumes differed between 
brain structures, and injury subgroups (GCS score, PTA, DAI) were asso- 
ciated with specific patterns of volume loss. Changes in cortical diffusion 
properties were detected in the more severe TBI subgroups in particular 
in patients with DAI grade 3, where diffusion properties increased in 
value and encompassed more cortical regions with time. Volume chang- 
es and diffusion changes were not associated. Cortical GM, hippocam- 
pus and lenticular nucleus volumes at 3 months were all very 
significant independent predictors of 12-month outcome even after 
adjusting for PTA duration, GCS scores and age. Thus TBI induced early 
and widespread structural brain volume loss which has significant im- 
pact on outcome in the chronic phase in moderate and severe TBI 
patients. 
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